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Abstract 

The value of the a spectroscopic factor (S a ) of the 6.356 MeV l/2 + state of 17 is 
believed to have significant astrophysical implications due to the importance of the 
13 C(a,n) 16 reaction as a possible source of neutron production for the s process. 
To further study this effect, an accurate measurement of the 13 C( 6 Li,d) 17 reaction 
at -E'lab = 60 MeV has been performed recently by Kubono et aJ., who found a new 
value for the spectroscopic factor of the 6.356 MeV l/2 + state of 17 based on a 
distorted wave Born approximation (DWBA) analysis of these data. This new value, 
S a ~ 0.011, is surprisingly much smaller than those used previously in astrophysical 
calculations (S a « 0.3 — 0.7) and thus poses a serious question as to the role of the 
13 C(a,n) 16 reaction as a source of neutron production. In this work we perform a 
detailed analysis of the same 13 C( 6 Li,d) 17 data within the DWBA as well as the 
coupled reaction channel (CRC) formalism. Our analysis yields an S a value of over 
an order of magnitude larger than that of Kubono et aJ. for the 6.356 MeV l/2 + 
state of 17 O. 

Key words: NUCLEAR REACTIONS, 13 C( 6 Li,d), £ lab = 60 MeV, DWBA and 
CRC analyses, deduced S a and reduced 7^ 
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1 Introduction 

The slow neutron capture, or s process is thought to be the production mech- 
anism for approximately half of all heavy elements in the universe [1,2]. The 
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asymptotic giant branch (AGB) phase of low and intermediate mass stars is 
thought to be the most likely astrophysical site for this process [2]. Although 
the neutron source for the s process has still not been well identified, the most 
recent studies seem to prefer the 13 C(o;,n) 16 reaction as the main neutron 
source in AGB stars at low temperatures [2,3]. 

Models of the s process depend critically on the neutron flux produced by the 
13 C(a,n) 16 reaction. However, the astrophysical S factor for this reaction, 
defined by: 

S(E) = a(E)Eexp(2nri) (1) 



where cr(E) is the cross section and rj the Coulomb parameter, has been de- 
termined experimentally down to a centre of mass energy of only 270 keV [4], 
whereas the reaction takes place predominantly at energies below this point. 

As the cross section for the 13 C(o;,n) 16 reaction decreases extremely rapidly 
as the incident a energy gets lower, direct measurement of the reaction rate at 
lower energies is very difficult. Extrapolations [4,5] suggest a rapid increase of 
the S factor as the incident a energy reduces to zero, which has been ascribed 
to the effect of resonances in 17 [5], the .F = 1/2" level at 5.94 MeV and 
the l/2 + level at 6.356 MeV, both very close to the a + 13 C threshold of 
6.359 MeV. However, due to the experimental uncertainties, the data are also 
consistent with a constant, horizontal extrapolated S factor [4]. 

In a recent paper, Kubono et ai. [1] suggested that a better way to determine 
the S factor at lower energies is via a direct a transfer measurement. The 
reaction chosen for this study was the 13 C( 6 Li,<i) 17 transfer at E\ a ^ = 60 MeV 
and angular distributions for transfers leading to the 0.0 MeV 5/2 + , 0.87 MeV 
1/2+ 3.055 MeV 1/2", 3.84 MeV 5/2", 4.55 MeV 3/2~ and 6.356 MeV 1/2+ 
states of 17 were accurately measured. The a spectroscopic factor S a of each 
state was determined from a finite-range distorted-wave Born approximation 
(DWBA) analysis [1] of these data. Using the value of S a obtained for the 
6.356 MeV 1/2+ state, its a width was determined and the cross section 
of the 13 C(a,n) 16 reaction through the tail of this sub-threshold resonance 
calculated. The astrophysical S factor was then calculated using Eq. (1). 

The DWBA analysis of Kubono et ai. gave a very small value of S a ~ 0.011 
for the 6.356 MeV 1/2+ state of 17 when normalised by the same scaling 
factor used to obtain S a = 0.25 for the 3.055 MeV l/2~ state. This leads 
to a very small contribution of the 1/2+ sub-threshold state to the reaction 
rate at low energies and thus to an astrophysical S factor that is essentially 
constant below £ c . m . = 300 keV. Since such a result poses serious questions 
about the mechanism of neutron production for the s process, a reanalysis of 
the Kubono et ai. data was undertaken. In the present work we present the 
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results of detailed DWBA and coupled reaction channels (CRC) analyses of 
these same data to determine whether S a for the 6.356 MeV l/2 + state of 
17 is actually so small. Our results show that the original DWBA analysis of 
Kubono et al. led to incorrect conclusions, and that their data are compatible 
with a much larger value of S a for the 6.356 MeV l/2 + state. The present 
analysis is consistent with previous extrapolations of the astrophysical S factor 
for the 13 C(o;,n) 16 reaction that indicate a rapid increase as the incident a 
energy approaches zero. 



2 The DWBA calculations 



Initially, we carried out a DWBA analysis as similar as possible to that of 
Kubono et al. [1], but using the code FRESCO [6]. We employed the same a 
+ 13 C binding potential and the same 6 Li + 13 C and d + 17 optical potentials 
as in the original analysis [1]. However, all calculations presented here were 
carried out including the full complex remnant term unless otherwise stated. 
The remnant term occurs in the residual interaction W in the expression for 
the transition amplitude in the DWBA. For the reaction A + a(— b + x) — > 
B(= A + x) + b the residual interaction is defined as follows for post form 
DWBA: 

Wp = V bB -U p = V bx + (V bA - Up) (2) 

and in the prior form: 

W a = V aA -U a = V xA + (V bA - U a ). (3) 



The remnant terms are the quantities in parentheses, U a and Up are the (com- 
plex) optical model potentials in the entrance and exit channels, respectively 
and V bx and V xA are the potentials binding the transferred particle (x) to the 
projectile and target cores, respectively. The quantity V bA is a (complex) op- 
tical potential operating between the projectile and target cores. See Satchler 
[7] for a full discussion of these terms. 

The use of the full remnant term enables good agreement between post and 
prior formulations of the DWBA to be obtained, which is not possible when no 
remnant term is included. We show results for the post formulation, although 
we did perform test calculations using the prior form to ensure that agreement 
was obtained. As a further check we compared calculations performed using 
FRESCO with no remnant term (using the post formulation) with calculations 
carried out with the code DWUCK5 [8]. For each transfer considered here 
excellent agreement between the two codes was obtained. 
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Table 1 

Values of N, L and S a used in the DWBA calculations. Values of S a are given for 
calculations with the full complex remnant term and calculations with no remnant 



term included. 


State 


r 


N 


L 


S a (full complex remnant) 


S a (no remnant) 


0.0 MeV 


5/2+ 


2 


3 


0.36 


0.15 


0.87 MeV 


1/2+ 


3 


1 


0.42 


0.18 


3.055 MeV 


1/2" 


4 





0.81 


0.32 


3.84 MeV 


5/2- 


3 


2 


0.64 


0.28 


4.55 MeV 


3/2- 


3 


2 


0.90 


0.38 


6.356 MeV 


1/2+ 


4 


1 


0.49 


0.24 



The entrance channel 6 Li + 13 C optical potential was set 1 of Kubono et al. 
[1] and we used their set 3 for the exit channel d + 17 O potential [9]. The a 
+ d binding potential was that of Kubo and Hirata [10], with the a particle 
assumed to be in a relative 2s state with respect to the d core. Throughout 
we adopt the convention that the number of radial nodes includes that at 
the origin but excludes that at infinity. The a + 13 C binding potential had 
a radius of 1.25 x (4 1/3 + 13 1/3 ) fm and diffuseness 0.65 fm, the depth being 
adjusted to give the correct binding energy for each 17 O state considered. 
The number of nodes iV and orbital angular momentum of the transferred a 
particle with respect to the 13 C core L were fixed by the oscillatory energy 
conservation relation 2(iV — 1) + L — J2t=i 2(n^ — 1) + k, where (rij, k) are the 
corresponding single-nucleon shell quantum numbers. The values of N and L 
used and the spectroscopic factors S a obtained in the analysis are given in 
Table 1. The spectroscopic factor for the a+d overlap was taken to be 1.0. 
Note that the S a values implicitly contain the C 2 term, where C is the isospin 
Clebsch-Gordan coefficient; however, in this case C — 1. As can be seen from 
Table 1 the contribution from the remnant term in the DWBA calculations is 
substantial and accounts for more than 60% of the S a strength for each 17 
state. The results of the DWBA calculations are compared with the data in 
Fig. 1. 

Figure 1 shows that we obtain reasonably good descriptions of the data, with 
the exception of that for transfers to the 3.84 MeV 5/2" and 4.55 MeV 3/2" 
states at forward angles. Unlike the original analysis of Kubono et al. [1], we 
obtain a reasonable description of the 0.0 MeV 5/2+ state; as the a particle 
is in an L = 3 state relative to the 13 C core the DWBA angular distribution 
should be relatively structureless, like the data. As we obtain such an angular 
distribution (both with FRESCO and DWUCK5), in contrast to the struc- 
tured angular distribution for this state shown in Ref. [1], we are forced to 
conclude that these authors have made an error in their calculations for this 
state, either through plotting the wrong curve or by using the wrong L value. 
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Fig. 1. DWBA calculations compared to the data. The dashed curves denote cal- 
culations including the full complex remnant term, while the dotted curves denote 
calculations with no remnant term. All calculations use the post form of the DWBA. 



It will be noted from Table 1 that if one adopts the procedure of Kubono et 
al. [1] of re-normalizing all the S a values with the scaling factor that gives a 
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value of S a = 0.25 for the 3.055 MeV l/2~ state, derived from a microscopic 
cluster calculation [11], one obtains from our DWBA results an S a value for 
the 6.356 MeV l/2 + state of 0.15 or 0.19, depending on whether the remnant 
term is included in the calculations or not. Either value is over an order of 
magnitude greater than that (S a = 0.011) found by Kubono et al. [1]. As we 
have used the same potentials in our DWBA calculations as those used in 
Ref. [1], we are unable to account for this discrepancy in our results for S a 
compared to theirs. Note that the use of the remnant term cannot account 
for the discrepancy (it is not clear from their paper whether Kubono et al. [1] 
included a remnant term in their calculations) as is apparent from Table 1. 
The only uncertainty is in the values used for the number of nodes N in the 
a + 13 C wave functions used by Kubono et al, as these are not given in their 
original paper. 

There is some ambiguity in the choice of N, depending on the structure as- 
sumed for the state in question. In deriving the N values for the 3.055 MeV 
1/2', 3.84 MeV 5/2 - and 4.55 MeV 3/2" states given in Table 1 we have as- 
sumed a 2p-lh configuration for these states. However, it has been suggested 
that these states are of 4p-3h character [12], and if one assumes this structure 
one obtains N values of 5, 4 and 4 for the 3.055 MeV 1/2", 3.84 MeV 5/2" 
and 4.55 MeV 3/2~ states, respectively. These larger N values lead to S a val- 
ues of 0.30, 0.25 and 0.36, respectively, which would yield a normalised S a 
value of 0.41 for the 6.356 MeV l/2 + state. In reality, these states are prob- 
ably a mixture of 2p-lh and 4p-3h configurations [11,13]; clearly, differences 
in the assumed structure of the 3.055 MeV l/2~ state cannot account for the 
discrepancy in normalised S a values between the current work and that of 
Kubono et al. [1]. 

For the 6.356 MeV l/2 + state itself, the value of N is more clearly defined; 
any physically reasonable choice of structure for this state results in a value 
of N = 4, the value we have adopted. Thus, we find that the smallest value of 
S a for this state consistent with the 13 C( 6 Li,d) 17 transfer data is 0.15 (after 
applying the normalisation procedure discussed above) and that it could be 
as high as 0.41, depending on the structure assumed for the 3.055 MeV l/2~ 
state, used to obtain the normalisation factor. It should be emphasised that 
there is no solid justification for the normalisation procedure used in Ref. [1] 
to obtain spectroscopic factors for all states based on the a-cluster structure 
of the 3.055 MeV 1/2 - state. Furthermore, this procedure is model-dependent 
as it relies strongly on the accuracy of the calculated S a of Furutani et al. [1 1] 
for the 3.055 MeV l/2~ state. 

We are thus forced to conclude that the original DWBA analysis of Kubono 
et al. [1] is flawed, and that the S a value obtained by them for the 6.356 
MeV l/2 + state is much too small. Our DWBA calculations using the same 
optical potential parameters as those used by Kubono et al. show that the 
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Table 2 

Parameters of the real DPP and volume WS imaginary potentials (5)- (6). The 
negative value for V indicates a repulsive potential. Radius parameters are given 
as: R x = r x x 13 1 / 3 . 

V (MeV) r v (fm) a v (fm) W (MeV) r w (fm) a w (fm) 



-18.0403 1.5079 0.6823 



37.2152 



1.3770 



1.0150 



13 C( 6 Li,ii) 17 transfer data measured by them are consistent with an (un- 
renormalised) S a ~ 0.49 for the 6.356 MeV l/2 + state of 17 0, well within the 
range used in astrophysical calculations (S a 0.3 — 0.7). 

We now discuss a subtle but quite important issue: the sensitivity of the 
DWBA calculation to the choice of optical potential for the entrance channel. 
This potential defines the door-way state of the transfer reaction and it should 
be determined as accurately as possible. The 6 Li+ 13 C elastic scattering was 
measured by Kubono et al. at forward angles only [1], which gives rise to 
some ambiguity in the choice of the optical potential. In Fig. 2 we plot the 
optical model fit to the 6 Li+ 13 C elastic scattering data using their parameter 
set 1. One can see that for large angles the predicted cross section exhibits 
unphysical behaviour compared to data for 6 Li + 12 C elastic scattering at the 
same 6 Li incident energy which extend to backward angles [14]. 

The 6 Li+ 12 C elastic scattering has been extensively studied in the past, and 
a very accurate systematics for the (energy dependent) optical potential has 
been established in a folding model analysis of these data over a wide range 
of incident energies [15]. An important result of this study is that the so- 
called dynamic polarisation potential (DPP) caused by 6 Li breakup has a 
strong repulsive contribution to the real potential. This DPP contribution 
can be rather well represented by a surface Woods-Saxon (WS) potential. We 
have, therefore, performed an additional optical model analysis of the 6 Li+ 13 C 
elastic scattering data at £i a b = 60 MeV using an entrance channel optical 
potential of the following form: 



U(R) = V Fold (R) + AV(R) + iW(R), 



where AV(R) = AVa v 
and W{R) = -W 



d 



dR 
1 + exp 



1 + exp 

R — R\y 

aw 



R-Ri 



ay 
-l 



(4) 
(5) 

(6) 



The DPP and volume WS imaginary potential parameters obtained (see Ta- 
ble 2) are quite close in shape to those determined from the Folding model 
analysis [15] of the 6 Li+ 12 C elastic scattering data at the same incident 6 Li 
energy but over a much wider angular range. 
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The Folding potential was calculated using the CDM3Y6 interaction (whose 
density-dependent parameters were fine-tuned to reproduce the bulk prop- 
erties of cold nuclear matter [16]) and the ground-state densities of 6 Li and 
13 C obtained in the independent-particle model by Satchler and Love [17] and 
Satchler [18], respectively. The predicted elastic scattering obtained using this 
potential is shown as the dotted curve in Fig. 2. Note that the volume integral 
of the real optical potential per interacting nucleon pair Jy (an important key 
to distinguish discrete potential families) is around -410 and -329 MeV fm 3 
for the set 1 potential taken from Ref. [1] and our Folding + DPP poten- 
tial, respectively. The empirical energy dependence of Jy is well-known from 
the global optical model analysis of light heavy- ion elastic data [19] , and one 
can find from this systematics that Jy should be around -330 MeV fm 3 for 
an energy of 10 MeV/nucleon (see Fig. 6.7 in Ref. [19]) which is in perfect 
agreement with our Folding + DPP calculation. We conclude, therefore, that 
the Folding + DPP potential should be more appropriate for the entrance 
channel optical potential. While both Kubono et al.'s set 1 and our Folding + 
DPP potential produce similar elastic scattering cross sections at the forward 
angles for which data are available, they differ considerably at larger angles, 
and a future measurement of the 6 Li+ 13 C elastic scattering at larger angles 
would be extremely helpful in determining a realistic entrance channel optical 
potential. 



We show the results of the DWBA calculations using this new 6 Li + 13 C 
optical potential in Fig. 3 and give the S a values obtained from this analysis 
in Table 3. As can be seen from Fig. 3, with the exception of the 3.055 MeV 
1/2 - state, the agreement between the DWBA calculations and the data is 
much better than that obtained using the set 1 potential of Kubono et al. 
This is particularly noticeable for the two l/2 + states. Table 3 also shows that 
the S a values obtained from this analysis are considerably smaller than those 
given in Table 1. The absolute S a value obtained for the 6.356 MeV l/2 + 
state is around 0.36 in this case. If we adopt the (controversial) normalisation 
procedure discussed above, we obtain S a values of 0.3 or 0.5 for the 6.356 
MeV l/2 + state, depending on whether a 2p-lh or 4p-3h structure is assumed 
for the 3.055 MeV l/2~ state. All these estimates of the S a value are over an 
order of magnitude larger than that (S a 0.011) determined in Ref. [1]. 



The influence of the entrance channel optical potential on the S a values ex- 
tracted from the DWBA analysis stresses again the importance of obtaining 
elastic scattering data over a sufficiently wide angular range for light heavy-ion 
systems such as that currently under consideration. 
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Fig. 2. Optical model calculations compared to the 6 Li+ 13 C elastic scattering data. 
The dashed curve denotes the predicted elastic scattering angular distribution using 
the optical model potential set 1 of Kubono et al. [1]. The dotted curve denotes the 
prediction of the present Folding + DPP potential. 



3 Coupled reaction channels calculations 



In order to test the possible influence of multi-step transfer paths on the S a 
values we also carried out a limited coupled reaction channels (CRC) analysis 
which included couplings between the ground state of 17 and its l/2f, 1/2]", 
5/2]", 3/2]" and 1/2^ excited states as well as the direct a transfer route. We 
did not consider multi-step paths proceeding via the excited states of 13 C as 
a physically meaningful analysis including these paths is not possible without 
some prior estimate of the necessary a strengths. 

We used the Folding + DPP optical potential described in the previous section 
in the entrance channel and parameter set 3 of Kubono et al. [1] in the exit 
channel. The 17 Coulomb coupling strengths were determined from the mea- 
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Fig. 3. DWBA calculations using the Folding + DPP optical potential for the 
6 Li+ 13 C system compared to the data. For the 3.055 MeV 1/2", 3.84 MeV 5/2" 
and 4.55 MeV 3/2" states the dashed curves denote calculations assuming a 2p-lh 
structure, while the solid curves denote calculations assuming a 4p-3h structure. All 
calculations use the post form of the DWBA. 
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Table 3 

Values of S a obtained from the DWBA and CRC calculations using the Folding + 
DPP optical potential for 6 Li+ 13 C system. 



State 


J 71 " 


N 


DWBA 


CRC (2p-lh) 


CRC (4p-3h) 


MeV 


5/2+ 


2 


0.15 


0.12 


0.15 


87 MeV 


1/2+ 


3 


0.17 


0.17 


0.17 


3 055 MeV 


l/2~ 


4 


0.30 


0.30 








5 


0.18 




0.18 


3.84 MeV 


5/2- 


3 


0.34 


0.34 








4 


0.19 




0.19 


4.55 MeV 


3/2- 


3 


0.48 


0.49 








4 


0.27 




0.27 


6.356 MeV 


1/2+ 


4 


0.36 


0.40 


0.40 



sured B(EX) values [20] and the nuclear coupling potentials were obtained 
by deforming the d+ 17 optical potential with nuclear deformation lengths 
obtained from the corresponding B(EX) values (assuming an 17 radius of 
1.2 x 17 1 / 3 fm). As the 17 O coupling strengths are weak, it was not found 
necessary to alter the <i+ 17 optical potential to compensate for the inclusion 
of these couplings. 

The results of the CRC analysis are shown in Fig. 4. The CRC description of 
the elastic scattering is almost identical to that shown by the dotted curve in 
Fig. 2. The S a values extracted from this CRC analysis are given in Table 3 in 
the two columns labelled "CRC(2p-lh)" and "CRC(4p-3h)", which denote the 
values obtained from the calculations assuming a 2p-lh and 4p-3h structure, 
respectively, for the 1/2]", 5/27/ and 3/27/ states. 

As Fig. 4 shows, with the exception of the transfer to the 5/2+ ground state, 
the shapes of the predicted angular distributions are unchanged by the inclu- 
sion of couplings between the ground state and the excited states of 17 0. The 
S a values are also either unchanged or, for the 5/2+ and 1/2+ states, very 
slightly altered by these extra couplings. The changes in the S a values for the 
5/2+ and 1/2+ states are smaller than those produced by, for example, the 
use of the Folding + DPP optical potential for the 6 Li+ 13 C system compared 
to Kubono et al.'s set 1 potential or by the inclusion of the full complex rem- 
nant term compared to a calculation with no remnant term. Therefore, we 
may conclude that couplings within 17 will have a negligible effect on the 
S a values extracted from an analysis of the 13 C( 6 Li,of) 17 transfer reaction. 
For the crucial 6.356 MeV 1/2+ state, an S a value around 0.4 is given by 
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Fig. 4. CRC calculations using the Folding + DPP optical potential for the 6 Li+ 13 C 
system compared to the data. The dashed curves denote the result of a calculation 
assuming a 2p-lh structure for the 3.055 MeV 1/2', 3.84 MeV 5/2' and 4.55 MeV 
3/2' states, while the solid curves denote the result of a calculation assuming a 
4p-3h structure for these states. 
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both DWBA and CRC calculations which is fully compatible with those used 
earlier (S a 0.3 — 0.7) in the s-process model calculations. 

However, multi-step effects can be important and lead to anomalous popula- 
tion of non a-cluster states in both the ( 6 Li,c?) and ( 7 Li,t) reactions, as was 
noted by Debevec et ai. [21] for the 12 C( 6 Li,<i) 16 reaction. They found that 
for 32 MeV incident 6 Li the 10.35 and 11.1 MeV 4+ states are populated 
with relative strengths of cr(ll.l)/cr(10.35) ~ 0.5 despite the fact that their 
a-particle widths are in the ratio r a (ll.l)/r a (10.35) xs 0.01. Kemper and 
Ophel [22] have shown that this anomalous population of the 11.1 MeV 4 + 
state is non-statistical in origin and must therefore be due to multi-step trans- 
fer paths. This problem is not so marked for the 12 C( 7 Li,t) 16 reaction and 
appears to become less important as the incident 6 Li or 7 Li energy increases. 
Nevertheless, it is an indication that multi-step paths can be important and 
must be considered carefully if one wishes to extract meaningful S a values 
from the analysis of a transfer reactions. 



4 Reduced a widths 



One may also derive the reduced a width, 7 2 , of a state from a DWBA analysis. 
The reduced a width is related to the spectroscopic factor and the a+ 13 C 
bound state radial wave function used in the DWBA calculations as follows 

ll = S a ^\u L (R)\ 2 (7) 



where /i a is the a+ 13 C reduced mass, R is the channel radius and ul(R) is the 
a+ 13 C bound state radial wave function, normalised such that / °° u 2 L (r)r 2 dr = 
1.0. We take R = 5.5 fm as our channel radius, the value used by Kubono et ai. 
[1] and Furutani et ai. [11]. Such a choice for the channel radius is reasonably 
compatible with the 'strong absorption' radius of the a+ 13 C system of around 
5.2-5.3 fm. The reduced a width is often scaled by the Wigner limit, defined 
as 

3h 2 

^(R) = ™ (8) 



to obtain the dimensionless reduced a width, 0^(R) [23]. It should be noted 
that Kubono et ai. [1] have made assumptions about the a+ 13 C bound state 
wave function that are equivalent to having S a = 9^, which is not always the 
case; we have not made this assumption here. 
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Table 4 

Q 2 a and 7^ values obtained from DWBA analyses using the 6 Li+ 13 C optical po- 
tentials of Kubono et al. [1] (Set 1) and the present work (Folding + DPP) at a 
channel radius R = 5.5 fm. Also shown are values for (9^,(5.5) from the microscopic 
cluster-model calculation of Furutani et al. [11]. 



State 


,r 


N 


el 


Set 1 
ll (keV) 


Folding 

01 


I + DPP 

II (keV) 


Ref. [11] 

el 


0.0 MeV 


5/2+ 


2 


0.209 


142 


0.087 


59.1 




0.87 MeV 


1/2+ 


3 


0.331 


224 


0.134 


90.8 




3.055 MeV 


1/2- 


4 


0.723 


490 


0.268 


182 


0.084 






5 


0.254 


172 


0.153 


104 




3.84 MeV 


5/2- 


3 


0.554 


375 


0.294 


199 


0.23 






4 


0.216 


146 


0.164 


111 




4.55 MeV 


3/2- 


3 


0.760 


515 


0.406 


275 


0.16 






4 


0.288 


195 


0.216 


146 




6.356 MeV 


1/2+ 


4 


0.272 


184 


0.200 


136 





In Table 4 we give values for 7^ and Q 2 a calculated using the S a values of 
Tables 1 and 3. Note that we do not use the normalised S a values in calculating 
■yl and 9l, as it is uncertain what the effect of the scaling process on the a+ 13 C 
bound state wave function should be. 

We have checked our procedures against those of Becchetti et al. [23] for the 
much more widely studied 12 C( 7 Li,t) 16 reaction and obtain good agreement 
with their values for 7^ for states in 16 using the same a+ 12 C wave functions 
and S a values. It should be emphasised that we obtain values of 7^ for the 
crucial 6.356 MeV l/2 + state considerably greater than the 7.4 keV value 
given in Ref. [1]. 

We also give in Table 4 6 2 a values for the 3.055 MeV 1/2", 3.84 MeV 5/2" and 
4.55 MeV 3/2~ states from the microscopic cluster-model calculation of Furu- 
tani et al. [11]. The calculations using Kubono et al.'s set 1 for the 6 Li + 13 C 
optical potential yield 8l values for the 3.055 MeV l/2~ state that are much 
smaller than the cluster model value, while for the 3.84 MeV 5/2~ and 4.55 
MeV 3/2~ states our Q\ values assuming a 4p-3h structure for these states are 
reasonably close to Furutani et aVs cluster model results. For the calculations 
using our Folding + DPP optical potential for the 6 Li+ 13 C system, the agree- 
ment is better and the cluster-model $1 values agree reasonably well with our 
values assuming a 4p-3h structure for these three states. This agreement with 
the calculation of Furutani et al. [11] shows that the d\ values that we have 
extracted from our DWBA calculations are at least plausible. Furthermore, 
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it supports the suggestion of Bethge et al. [12] that these states (the 3.055 
MeV 1/2- , 3.84 MeV 5/2 - and 4.55 MeV 3/2 - ) are of predominantly 4p-3h 
character. 



5 Conclusions 

A recent DWBA analysis of new 60 MeV 13 C( 6 Li,<i) 17 transfer data [1] found 
a very small a-spectroscopic factor (S a ~ 0.011) for the 6.356 MeV l/2 + state. 
The authors of Ref. [1] concluded that there is no large enhancement of the 
13 C(a,n) 16 reaction rate at energies of astrophysical interest due to this sub- 
threshold state based on the small reduced a width they obtained, contrary 
to previous suggestions [4,5]. Since such a result poses serious questions about 
the 13 C(o;,n) 16 reaction as a neutron source for the s process, we have per- 
formed DWBA and CRC analyses of their data using a realistic choice for 
the entrance channel optical potential. We concluded that the transfer data 
are consistent with a much larger S a value for the 6.356 MeV l/2 + state and 
consequently a much larger 7^ value. Thus, the transfer data are compatible 
with a large contribution from this state to the 13 C(a,n) 16 reaction rate at 
very low energies. 

We have also shown that uncertainties in the assumed a+ 13 C structure of the 
17 states can have significant effects on the extracted spectroscopic factors 
and reduced a widths. While this problem does not directly affect the 6.356 
MeV l/2 + state, one needs to be aware of it. As relative S a factors are rea- 
sonably well determined one may use the normalisation procedure employed 
by Kubono et al. to scale the DWBA results to a calculated value for a state 
that is believed to have significant a clustering structure in order to obtain 
absolute S a values. However, such a method is highly questionable as the frag- 
mentation of the a-cluster states over the low-energy region of 17 excitations 
is unknown. Moreover, there are difficulties in this method as it relies on the 
accuracy of the structure calculation used for the normalisation, plus it is not 
clear what effect this process should have on the a + core wave function if 
one wishes to extract 7^ values. 

We also noted the effect of the entrance channel optical potential on the S a val- 
ues extracted from the DWBA analysis, which was considerable. This demon- 
strates the need to have a realistic choice for the optical potential in the 
entrance channel. It is, therefore, necessary to obtain elastic scattering data 
over a large angular range (possibly much larger than the angular range of the 
transfer data themselves) if one wishes to use this type of indirect approach 
to obtain reliable estimates of astrophysical reaction rates. 

A limited coupled reaction channels study found that couplings between the 
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ground and excited states of 17 had a negligible effect on the S a values. Other 
multi-step transfer paths, in particular, those proceeding via the excited states 
of 13 C, were not investigated due to the need for some prior estimate of the 
necessary a strengths. Without such estimates one merely introduces extra 
adjustable parameters that need to be determined from the same data set. 

In summary, we have highlighted some of the pitfalls of DWBA analyses of a 
transfer reactions and their use to extract spectroscopic factors. We have not 
considered in detail the more fundamental question of whether such transfer 
reactions may be adequately described by a simple direct a transfer, as is 
assumed in the DWBA. If multi-step paths contribute significantly to the a 
transfer strength the S a values extracted from such a DWBA analysis will not 
be reliable, regardless of how carefully it is carried out. However, an extensive 
investigation of these multi-step transfer paths still cannot be carried out in 
a meaningful way for the 13 C( 6 Li,cf) 17 reaction due to the lack of prior esti- 
mates of the required a strengths for transfers proceeding via the excited states 
of 13 C. Thus, while the indirect approach to astrophysical reaction rates can 
be a valuable tool one needs to be aware of possible uncontrolled parameters 
contained in the analysis. 
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